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We present analytic and numerical results for the thermoelectric efTect in unconventional super- 
conductors with a dilute random distribution of impurities, each scattering isotropically but with a 
phase shift intermediate between the Born and unitary limits. The thermoelectric response function 
has a linear temperature dependence at low temperatures, with a slope that depends on the impurity 
concentration and phase shift. Although the thermoelectric effect vanishes identically in the strict 
Born and unitary limits, even a small deviation of the phase shift from these limits leads to a large 
response, especially in clean systems. We also discuss possibilities of measuring counter-flowing 
supercurrents in a SQUID-setup. The non-quantized thermoelectrically induced flux can easily be 
of the order of a percent of the flux quantum in clean systems at ''ife temperatures. 

PACS numbers: 74.25.Fy, 74.72.-h 



I. INTRODUCTION 

Electronic charge and heat transport measurements 
can give important information about the micro- 
scopic properties of high-Tc and other unconventional 

SUperCOnductorS|li^i3,4,5,6,7,8,9,10,ll,12,13,14^15,16,17,18,19,20 

Of particular interest is the low-temperature regime, 
where electronic properties are believed to be controlled 
by elastic impurity scattering. This is the case, because 
in an unconventional superconductor with an order 
parameter changing sign at some points on the Fermi 
surface (for example the d-wave order parameter) , elastic 
scattering by even a small concentration of impurities 
leads to pair breaking, and the formation of an impurity 
band of width 7 around the Fermi level. These low- lying 
excitations, are responsible for a universal limit at 
low-temperatures T ^ 7, where the charge and heat 
conductances become independent of the scattering 
properties of the impurities. This low-energy behavior 
was first predicted by hee— for the zero-frequency 
charge conductance, and was later extended to include 
the heat conductance by Graf et al.m^ 

The universal character of the T ^ hea t conduc- 
tance was seen experimentally in Rcfs. [^iJlaallyllll • A 
series of recent related experiments on the temperature 
dependence of the heat conductance approaching the uni- 
versal limit ^ microwave conductivity^ and STM spec- 
troscopy of quasiparticle impurity states;^ have all in- 
dicated that impurity scattering in the cuprates might 
not be in the strict unitary or Born limits. Limiting 
the discussion to isotropic scattering, this means that 
the impurity potential uq is of intermediate size, and 
the corresponding s-wave scattering phase shift So = 
arctan(7rA/'/uo) has some intermediate value < So < 
7t/2. Here is A// the density of states at the Fermi level in 
the normal state. If this is the case, electron-hole symme- 
try is explicitly broken near each impurity and the global 
electron-hole symmetry is broken for a homogeneous di- 
lute distribution of such impuritiesiSSiSiSiS^ This has 
direct consequences for the heat and charge transport 



coefficients but leads also to large thermoelectric ef- 
fects. A fact that was previously noted in connection 
to the heavy fermion systems by Arfi et alm^t^ In the 
present paper we report an extensive analysis of electron- 
hole symmetry breaking by elastic impurity scattering, 
and its effect on the transport coefficients in unconven- 
tional superconductors, with an emphasis on the d-wave 
cuprates at low temperatures. 

Previous work on thermoelectric effects in the heavy 
fermion superconductor322i2liS& were limited to situations 
in which the energy-dependent broadening of the quasi- 
particle states can be neglected, i.e. when the imaginary 
part of the impurity self-energy satisfies 5S|^(e) ^ e. 
This approximation is expected to be good at high tem- 
peratures, T < Tc, but fails at lower temperatures. In 
fact, in the temperature regime T < 7 of high current 
interest, where effects of universality become of impor- 
tance, the impurity renormalization is the dominant en- 
ergy scale and the problem has to be considered anew. 

Our results can be summarized as follows. (1) For 
intermediate phase shifts, the thermoelectric response 
function ri{T) is in general large, which leads to counter- 
flowing supercurrents detectable e.g. as a thermally in- 
duced magnetic flux in a ring setupi^LSSiSi^fli^ (2) At 
low temperatures T < j, ?7(T) scales linearly with tem- 
perature, with a non-universal slope that grows large in 
clean systems. 



II. ELECTRON-HOLE SYMMETRY BREAKING 
AND GIANT THERMOELECTRIC EFFECTS 



To compute the response to a thermal gradient VT and 
an electric field E{uj), we need to evaluate the charge 
current Sje as well as the heat current Sj^. In the lin- 
ear response, the observable response functions such as 
the heat and charge conductivities, can conveniently be 
expressed in terms of response functions ha/s that are 



2 



defined as 

('i)^(h hA( i V (1) 

^ ) \ L21 L22 j \ V(^) ) ^ ' 

This particular choice of forces ensures that the Onsager 
relations L12 = L21 are fulfilled (see e.g. Ref. [s^). 

Contrary to the two diagonal terms Ln and L22 in 
Eq. 1^, the thermoelectric coefficients L12 = L21 require 
an electron-hole asymmetry around the Fermi energy in 
order to be non-zero. In the normal state, to quasiclassi- 
cal accuracy, i.e. to leading order in small parameters 
s — {T/Tf,l/pf£^o,...}, there is no such electron-hole 
asymmetry in the theory and the thermoelectric coef- 
ficients vanish, = = 0. Here is ^0 = Vf/Tc the 
superconducting coherence length, and T/, p/, and Vf 
are the Fermi temperature, momentum and velocity, re- 
spectively. This result holds also in the superconducting 
phase when the order parameter has the conventional 
s-wave symmetry^ However, for unconventional super- 
conductors this is not the case: impurity scattering is 
pair-breaking and potential scattering off an impurity site 
induces a non- vanishing electron-hole asym metry already 
to leading order in the small parameters s |22l23l24l25l | . 
To illustrate this, we consider a superconductor having 
an order parameter with a vanishing Fermi surface aver- 
age, J dpfA(pf) = 0, in which case the equilibrium-state 
impurity fj^p-matrix in Nambu (electron-hole) space has 
the form 

£fl ^ sin(5o cosSoi + sinSo J dp fg^{pf,e)T3 
™P ttTV/ cos^^o - sin^i^oi/ dpfg^{pf, e)]^ ' 

where the diagonal component of the equilibrium qua- 
siclassical Green's function is = — e^/J7^ [fl^ = 
•y/|A(p/)p — (e-'^)^]. The electron- hole asymmetric scat- 
tering becomes explicit when we examine the electron 
(11) and hole (22) parts: 

^"(^'^ ^ cos6oTsmSoJdpfg^iPf,ey 

In particular, multiple scattering off the impurity leads to 
resonances at e = —ercs{6o) for electron-like excitations 
and at e = +£^03(60) for hole- like excitations. These im- 
purity resonances become virtual bound states localized 
at the impurity in the strong scattering limit, Sq — > 7r/2 
[33^ . This also implicates that electrons and holes have 
different energy dependent scattering life times at inter- 
mediate phase shifts 

We assume that the impurities are randomly dis- 
tributed with an average small concentration n^, that 
is absorbed in the normal state elastic scattering rate 
TQ^ = 2r~ 2ni sin^ So/TrMf- Within the usual impurity 
averaging technique^ the impurity self energy has the 
form 

tl^ie) = S]^(6)i + (e)f3 = ndtpi^). (4) 
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FIG. 1: The imaginary part of Efi shown for different nor- 
mal state scattering phase shifts So for a d-wave supercon- 
ductor with A{pf) = Ao cos(2(^p'j, ). There is a substantial 
asymmetry in Q'Efi(e) around e = 0. The dashed vertical 
lines indicate the single-impurity r6SOI13,IlC6 6IlGrgl6S 6rcs 
for an electron-like quasiparticle at the corresponding phase 
shifts. The resonance energy of the single impurity self energy 
is carried over to the self-consistently computed self energy 
describing the randomly distributed impurities, but slightly 
shifted away from the Fermi level; e(5o) ^ eres(<5o). The 
hole component, E22, is related to the electron component 
asE«(6) = Eo«-Ef = [Efi(-6)]*. 



The term proportional to the third Pauli matrix in 
electron- hole space, fa, gives the impurity renormaliza- 
tion of the energy, = e — T,f . The unit term, So^(e), 
drops out of the equations for the equilibrium state, but 
in fact enters the transport equations describing the non- 
equilibrium state. S^(e) explicitly breaks electron-hole 
symmetry, which leads to the thermoelectric effect we 
study below. In general we can not expect the unit term 
to either vanish or to be small, although it vanishes in the 
strict Born and unitary limits. The induced asymmetry 
is shown for the electron part E{\ = Eg' -I- in Fig. ^ 
The response functions hap, calculated to lowest or- 
der in the small parameters s, are conveniently computed 
through the quasiclassical Keldysh propagator Sg^ , eval- 
uated to linear order in the forces. Amazingly, 6g^ has a 
closed form in which the self-consistently computed equi- 
librium Green's function g^ , order parameter A^, and 
impurity self-energies S^-^p serve as input ^S*^ A compli- 
cation is that Sg^ depends on linear corrections to the 
self-energies 5a^ and 5a^ , which need to be computed 
self-consistently with 6g^ . This procedure is equivalent 
to take into account vertex corrections in the Kubo for- 
malism. In our model, we assume isotropic scattering, in 
which case the vertex corrections vanish. This is a well 
known fact, and is ultimately due to the anti-symmetric 
form of the forces Vf ■ E and Vf ■ VT. For cases of 
anisotropic impurity scattering (p-wave etc.), these cor- 
rections must be taken into account, as was done for the 
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diagonal response functions Ln and L22 by Durst and 
Leei^ In the following we neglect anisotropic scattering 
and compute all response functions L^p. 

The expression for Sg^ given in the appendix of 
Ref. 01 is general enough to serve as the starting point 



for our calculation. In fact, the calculation of the zero- 
frequency limit of the off-diagonal response functions 
L12 = L21 follows the same logical steps as the calcu- 
lation of Lii and L22 in Ref. [13 [their Eqs. (29) and 
(30)], and we only give the final result: 



T u 

^12 



de e sech 



2T 



AA(p>,e)3I]^(e) 



[5Rf7«(p;;6)]^-[3S«(e) 



(5) 



r 



/VL^[pf\ e)] is the density of 



where Af{pf , e) 

states in the superconducting phase. L12 is directly pro- 
portional to the imaginary part of the unit term 3 Sg- 
of the equilibrium impurity self energy, which is an odd 
function of energy. Eq. ^ is the proper generalization to 
arbitrary low temperatures, including impurity scatter- 
ing renormalization effects, of the results in Refs. p3l23 |. 



Conductivities 



Once we know L^^, 



we may compute the conductivi- 
ties. In the following we consider transport in the cuprate 
planes, along one of the anti-nodes of the c?-wave order 
parameter. We can then drop the superscripts of L^^- 
The charge and heat conductances are given by 



aiT) 
KiT) 



UijT) 
T ' 

L22(T) 
y2 



(6) 



In the presence of a temperature gradient, the non- 
vanishing thermoelectric coefficient L12 induces a charge 
current Sj^ = —rj^T, where we define 



viT) 



Li2(r) 

rp2 ■ 



(7) 



The appcarcncc of a bulk charge current, 6je, leads to 
a counter-fiowing supcrcurrent, which we discuss in the 
next section. 

In Fig. 121 we present the thermoelectric response func- 
tion 77 as a function of temperature for several phase shifts 
for a fixed rather short mean free path £ = vjTq ~ 5^0j 
where ^0 = / /Tc is the superconducting coherence 
length. In our model, the physical transition temper- 
ature Tc is then suppressed by about 25% compared 
to the clean limit transition temperature Tcq, in ac- 
cordance with the Abrikosov-Gorkov form for the Tc- 
suppression as function of F = (1/2to) by elastic impu- 
rity scatteringi^ vi'T) is sizable over a wide temperature 
range from zero to but vanishes in the T ^ Tc limit, 
where the electron-hole asymmetry is of order and 
neglected in our quasiclassical theory. The thermoelec- 
tric effect is sensitive to the microscopic superconducting 



properties, such as the order parameter size and its sym- 
metry, through the coherence factors entering Eq. |SJ)- It 
is also sensitive to the nature of the impurity scattering: 
first to the normal state mean free path, but also to the 
phase shift (5o. In the strict Born {5q — > 0) and unitary 
((5o 7i'/2) limits, r]{T) vanishes since the electron- hole 
asymmetry vanishes in those limits, 9S^(e) = 0. How- 
ever, for all other values of Jq, rjiT) is large, and has a 
maximum of the order of 0.1(7Ar/e near T ~ 0.5Tc. Here 
is (Jn = e'^Affv'j / {2T) the Drude conductivity in the nor- 
mal state. In fact, even a small deviation of the phase 
shift from e.g. the unitary limit leads to dramatic changes 
in the thermoelectric response, see below. At low tem- 
perature, we find 77 cx T, (see left half of Fig. The 
slope rj/T contains detailed information about material 
parameters such as the scattering phase shift, as we will 
discuss in more detail in Section IlIII 

We should mention that we have neglected inelastic 
scattering in our calculations. In e.g. Ref. 0| inelastic 
scattering by anti-ferromagnetic spin fluctuations were 
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FIG. 2: The thermoelectric coefficient 77 scaled by its low- 
T dependence r) ^ T for the scattering rate F = O.lTc and 
several normal state scattering cross sections ai = s\v? 5q. 
The inset shows the unsealed function. In the universal limit 
{T <^ 7), rj/T approaches a non-universal constant given in 
Eq. lO. 
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FIG. 3: A sketch showing the major features needed for the 
detection of thermoelectricaUy induced magnetic fluxes. The 
loop consists of two arms, one being an s-wave superconductor 
for which ry is orders of magnitude smaUer than in the d-wave 
arm. The two superconductors are connected by Josephson 
junctions and the d-wave arm is connected to two bulky re- 
gions were the temperature is regulated. 



included in the calculation of the thermal conductivity, 
and it was shown to give rise to the charac teristic increase 
in k{T) seen in experiments (e.g. jllll2l |'l just below Tc- 
Thus, we expect corrections to our results in Fig.|21in the 
high temperature region T <Tc. In the important low-T 
region, T ~ 7, inelastic scattering is not of importance, 
since the corresponding self energy scales as and is 
small compared to the self energy from elastic impurity 
scattering. 



B. ThermoelectricaUy generated magnetic flux 

In the normal state, the thermoelectric effect leads to 
the appearance of a voltage (given by the thermopower). 
However, in the superconducting state, a steady state 
voltage, and an associated electric field, is short circuited 
by the appearance of a supercurrent and an associated 
phase gradient (2L22i2S 



{e/m)nsPs = -Sje 



(8) 



where is the superfluid density and ps = ^(^X^ if ■'^) 
is the superfluid momentum. Thus, the total charge cur- 
rent is zero, j^o***' = Sje + is ~ 0, but the phase gradient 
can be dete cted in a flux measurement as was done in 
Ref. Isfiaij . These early experiments were carried out 
with low- Tc s-wave superconductors, for which thermo- 
electric effects are smaller than in the unconventional su- 
perconductors we are considering, for two reasons. First, 
they are down by the electron-hole asymmetry factor s; 
second, they are exponentially suppressed at low tem- 
peratures by the energy gap around the Fermi surface. 



FIG. 4: The magnetic flux induced by a temperature gradi- 
ent in a ring with a thermoelectric response given in Fig. |5| 
It is assumed that a temperature difference 5T — Ti — T2 = 
O.OOSTc is maintained over the loop in Fig. |5] Current ex- 
perimental precision for flux measurements is 10~'^$o /VHz. 
Part (b) shows the divergence of the penetration depth in the 
d-wave superconductor at temperatures near Tc. 



Consequently, the experiments were mainly done in the 
temperature region near Tc, where the temperature de- 
pendence of Ug plays an important and somewhat para- 
sitic role when the goal is to measure 77 (T) (for a recent 
discussion of the experimental situation see Ref. H^)- 

To measure the flux generated by the counter-flowing 
supercurrent, we propose a hybrid SQUID-hke setup, as 
indicated in Fig.|31 with one arm of a low- Tc material and 
the other of the cuprate material of interest. The two 
arms are in electrical contact via two Josephson junc- 
tions. Thus, by construction, the generated flux will 
mainly originate from the cuprate arm and we can pre- 
dict a non-quantized flux of size 



<&(T) 



n2'K 



\{T f e7j{T) ST 

CTN F 



\2 
-^0 



(9) 



Here is $0 = co/2e the flux quantum, cq the ve- 
locity of light, ST = Ti — T2 the temperature bias, 

Ao = y^CQ / (iTTc'^Aff Uy- ) the zero-temperature penetration 

depth in the clean system, and A(r) = c^m/ {'i-Ke^ng) 
the temperature dependent penetration depth in the 
dirty system. We assume that the equilibrium flux is 
zero and put the integer n = hereafter. The super- 
fluid momentum is proportional to the local temperature 
gradient VT, which leads to that the flux $ is propor- 
tional to the temperature bias ST, since $ is related to 
a contour integral of ps around the loop (the evaluation 
of (f> is a standard calculation that can be found in e.g. 
Ref. lU ) • Note that if both arms were of the same ma- 
terial, the thermoelectric response in the two arms would 
give counter-flow in opposite directions around the loop 
and cancel. It is therefore necessary to have different 
responses in the two arms, although it does not matter 
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FIG. 5: The thermoelectrically induced flux as a function of temperature for strong scattering impurity potentials, (5o = 75° (pi « 
0.93) and 5o = 85° (ai « 0.99), at different scattering rates (the corresponding mean-free-paths are: 5,50, and 500 ^o)- In the 
two left panels we show the actual flux generated as a function of temperature, setting liT — 0.005Tc. At high temperatures, 
T <Tc, the flux is a sizable fraction of $o and a larger ST may well give a flux > $o- At low temperatures T < O.lTc the flux 
may still be of order 0.1% — 1% $o if the material is clean, £ > 20 ^o, as indicated in the insets. In the two panels to the right 
we show the factor (\{T)/\o)'^er]{T)/aN ■ This factor depends strongly on the scattering phase shift So (see also Fig.|5|where 
erj{T)/aN is displayed for a larger variety of scattering phase shifts) but does not have a particularly strong dependence on F 
apart from close to Tc where the temperature dependence of X{T) dominates. 



how that is achieved. Our setup is Hmited to low temper- 
atures by the low transition temperature ^ 10 — 15 K of 
the s-wave superconductor. To probe the high temper- 
ature regime of Eq. Q, some other type of asymmetry 
in the thermoelectric response of the two arms must be 
accomplished. Disregarding these complications, we now 
investigate Eq. (O for all temperatures below Tc of the 
d-wave arm. 

In Fig. 0] we give an estimate of the induced flux cor- 
responding to the thermoelectric response in Fig. [21 as- 
suming a temperature difference of O.OOSTc sa 0.5K main- 
tained across the loop. When the scattering phase shift is 
intermediate between the Born and unitary limits, fluxes 
of order of 0.1% — 1% of the flux quantum is generated 
at ^He temperatures ~ A.2K (~ 0.05Tc). At higher tem- 
peratures the temperature difference could be allowed to 
be larger and we can predict fluxes on the order of 10% 
of $o- For cleaner samples, the induced flux can be even 
larger, of the order of a few percent of $o at T = 0.05 T^ 
see Fig. |31 

For high temperatures, the temperature dependence 
of the flux is heavily influenced by the temperature 
dependence of the penetration depth, see Fig. 0] and 
Fig- El Thus, to extract the T-dependence of i]{T) from 
<i>(T), the temperature dependence of A(T) should be di- 
vided out. However, in the important low-temperature 
regime (to be discussed further in the following sec- 
tion), $ is directly proportional to ri(T) cx T, since the 
penetration depth is limited by impurity scatteringjSi^ 
A(T) — A(0) cx T^, and its T-depcndence can be neglected. 



We note that the minimum temperature at which mea- 
surements can be performed will be set by a combination 
of two factors: first, the smallest temperature bias that 
can be applied (ST ^ T is needed in order to have a uni- 
form thermoelectric response function 77 (T) in the sam- 
ple); second, the flux measurement sensitivity, since the 
flux scales as $ - rj{T)ST ~ TST. Thus, at low tem- 
peratures, there is a trade off between having a small 6T 
and at the same time have a measurable flux. 

Note that the induced superfluid momentum ps is 
small. We estimate 



Tr 



L $0' 



(10) 



where L is the distance between the two reservoirs (i.e. 
we have set VT = 5T/L). Thus, we do not need to take 
feedback effects from the superfiow via the Doppler shifts 
into account. 



III. LOW TEMPERATURES 

In the low-T limit, the response functions can be 
studied analytically through a systematic Sommerfeld ex- 
pansion in the small parameter T/7. The magnitude of 
7 depends on the normal state scattering rate as well 
as the phase shift, and we have 7 ^ ^J'^FKQ^J2 in the 
unitary limit while 7 ~ 4Ao exp[— 7rAo/2r] in the Born 
limit. For intermediate phase shifts 7(r, 5q) can be found 
numerically as we discuss below. Only in the limit T ^ 7 
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FIG. 6: (a)-(c) Phase-shift dependence of the impurity band width 7 and the parameters a and c normalized by "i/Tc, see 
Eqs. lAlll - IIASII . (d) The thermoelectric coefficient 77/r at low temperatures T ^ 7. The dashed (red) curve contains the low-T 
constants in Fig.|5] We have used A(p/) = Ao cos 2(j!)p-^ throughout. 



is the Sommerfeld expansion useful. Thus, wc consider 
low energies and write 



~R 



^ ae + + be^) + de" + 0[t% 



= ce + he" 



(11) 



with real constants a, 6, c, d, h, and 7. These constants 
are determined self-consistcntly by the low-energy equa- 
tions presented in the appendix. Inserting this ansatz 
into Eq. ||SJ), keeping terms to leading order in e, we find 
for the thermoelectric response 
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were ^ = (l/Ao)|dA(0)/o?(/)|0„^^^ is the opening rate of 
the gap at the nodes (related to the so called gap velocity, 
see Ref. [S^). The coefficient of the correction to the 
T — > asymptotic is given by 



ai2 



2c' 



3/172 



357. 



(13) 



Including phase shifts away from the unitary and the 
Born limits will also change the rate at which the charge 
and heat conductivities approach their universal values. 
Wc find the leading order dependence of the charge con- 
ductivity to be 



CT(r< 7,W"> 0) 



ttmAq 

and of the heat conductivity to be 



2 2 T^2 ^ 
1 + — ^i^^l, (14) 
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where the coefficients of the T^-terms contain direct 
information about the impurity induced particle-hole 



asymmetry through the parameters a and c, 



^22 



2c' 



(16) 



In the zero temperature limit T ^ 0, both the charge and 
heat conductivities approach universal values (indepen- 
dent of the properties of the impurities) while rj/T has a 
non-universal T —^ limit set by the ratio c/7. All func- 
tions, including 77, have non-universal low-T corrections 
to their T = values and are sensitive to the phase shift 
(5o and the scattering rate F. We note that the low-T ex- 
pressions given in Ref. |l8j for a and k/T were implicitly 
restricted to the unitary and Born limits. Consequently, 
the parameter c is absent in their Eqs. (56)- (57). 

We can analyze these results further by solving the 
self-consistency equations for the impurity self energy to 
lowest order in e, and thereby determine the parame- 
ters 7, a, c, etc. The results are given in the appendix, 
Eqs. (|A1|) - (|A6|) . and the numerical solution is presented 
in Figs. l6l7l The impurity band width 7 is exponentially 
small for phase shifts far from the unitary limit. As a 
consequence, the universal limit is reached at an expo- 
nentially small temperature. This is also confirmed in 
the left half of Fig. |21 where we present the thermoelec- 
tric response function on a logarithmic scale. Clearly, the 
crossover to the low-T power law is severely pushed down 
in T for small phase shifts. This suppression becomes ex- 
tremely fast when the system is clean. In fact, for ultra- 
clean samples, with F ^ lO^'^Tc or lower (corresponding 
to normal state mean free paths £ ~ 500^o or longer), a 
deviation of the phase shift away from 7r/2 by only a few 
percent will reduce 7 by several orders of magnitude, see 
Fig- da). At the same time, the ratio c/7 is increased 
humongously [Fig. El^c)], t]/T grows large [Fig.[7|^c)] and 
the thermoelectric coefficient itself is sizable on the scale 
of ajv/e [Fig. Eld)]. 
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FIG. 7: (a) The impurity band width 7 is extremely sensitive to phase shift variations in clean systems, (b) The ratio of the 
parameters c and a. (c) The low-T asymptotic of the thermoelectric coefficient grows large even for small deviations of the 
phase shift from the unitary limit, (d) The temperature dependence of the thermoelectric response function for a phase shift 
So ~ 86.4° ((Ti = sin^ Sq = 0.996). Note the data collapse when rj is scaled by aN/e, because 77 grows as l/F in clean systems. 



The thermoelectric response scales as l/F, which is 
also clear in Fig.jS] Thus, if wc change the normalization 
in Fig.EIc)-(d) from the normal state Drudc conductivity 
(Tat (which effectively cancels this scaling), to the univer- 
sal limit charge conductance a{T = 0) = (Too in Eq. (|14|l . 
we get the scale Tceri/[Taoo) ~ 100 - 1000 in Fig. EJc), 
and CTy/croo - 1 - 10 in Fig.|7|;d), for F - 0.01 - O.OOlTc- 
Thus, the thermoelectric effect grows large in clean sys- 
tems with phase shifts close to the unitary limit. 

Note the equal importance of the parameters c and 
a in a region near the unitary limit, see Figs. EJb)- 
(c) and Fig. Ob). The ratio approximately behaves as 
c/a ^ (7/F) sin(5o cos Jo- Thus, unless 7 ^ F (which 
holds far from the unitary limit deep inside the region 
where 7 is exponentially suppressed) it is in general not 
allowed to neglect the unit term when conductivities are 
computed at intermediate phase shifts, although that has 
been common in the literature (e.g. in [112011361 ). 

Finally, by comparing Figs. |H| and |21 we see that the 
true universal limit is in fact reached at very small tem- 
peratures, T < 10~^7. This is due to the largeness of the 
parameters a/7 and c/7 [c.f. Fig.|Bfb)-(c)] and higher or- 
der parameters 6, d, and ft,, for clean systems with phase 
shifts near but not strictly equal to 7r/2. Thus, we need 
T < il/o.) ^ 7 in order for the Sommerfcld expansion 
to work well. Or, strictly speaking, we have an effec- 
tive Sommerfcld expansion parameter T/{'-f/a), instead 
of r/7. 

IV. SUMMARY AND CONCLUSIONS 

In this paper we have focused on phase shifts close 
to 7r/2, because that appears to be the relevant re- 
gion experimentally. Perhaps the most clear signature 
of this is the impurity induced resonances seen in STM 



experiments^ The resonance energy is near the Fermi 
level, but not at the Fermi level, thus signaling a phase 
shift close to 7r/2. Also, if the low temperature heat con- 
ductivities found in cxperimentsSiiiSiSiifliiiii^ are indeed 
universal, the phase shift needs to be close to tt/2. Oth- 
erwise, 7 is exponentially suppressed 7 F and it would 
be hard to argue that the limit T ^ 7 was reached. The 
limit i5o = arctan(7rA//Mo) — *■ 7r/2 can in real life of course 
only be approached asymptotically. It is therefore likely 
that thermoelectric effects are of importance in the high- 
Tc, as well as other unconventional superconductors (uni- 
versal heat conductance was also found in Sr2Ru04 |^). 
The effect grows large in clean systems. The ideal exper- 
iment would then be to measure the thermoelectrically 
induced flux, together with the thermal conductivity, at 
low temperatures for samples with different levels of dis- 
order, e.g through controlled Zn doping. The slope of the 
thermal conductivity is related to the gap size Aq and 
its slope at the node /i, through (7r^/3)(2A//Uj)/(7r^Ao), 
which can be compared with values of Aq and ^ from 
other experimental probes (e.q. photoemission^) . The 
ratio rjl{eK) = c/7 contains information about F and 
(5o through Eqs. HAip - HA3|l . with Aq and as input pa- 
rameters. We note that the exact angle dependence of 
the gap function is not particularly important and it is 
enough to know Aq and /i. If the normal state scattering 
rate F is known independently, the scattering phase shift 
is uniquely determined by 77/ (bk). 

When the energy scale 7 is exponentially suppressed 
in clean systems, the temperature region in which the 
term in the expansions in Eq. H12f) - (|15|1 is of impor- 
tance, is also suppressed. Thus, it is not necessarily ideal 
to study super-clean samples. Another issue with super- 
clean samples is the validity of the homogeneous scatter- 
ing model with only s-wave scattering, in which a dilute 
concentration of point impurities are assumed to be ran- 
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domly distributed in the sample. For sufficiently clean 
samples this model has to break down and should be re- 
placed by more realistic models of extended impurities 
(suggestions can be found in e.g. Refs. 

In conclusion, we have computed the thermoelectric 
effect in unconventional superconductors with impuri- 
ties scattering in neither the Born nor the unitary lim- 
its. The thermoelectric effect is an interesting unex- 
plored avenue for the investigation of microscopic prop- 
erties of high-Tc as well as other unconventional super- 
conductors, in particular at low temperatures. Of spe- 
cial interest is to extract material parameters such as 
the gap size Aq, the slope of the gap at the gap node 
fjL = (l/Ao)|(iA(0)/(i(/)|0^^^^, the impurity band width 
7, and the impurity scattering rate F and phase shift 
5o. This can be accomplished by measuring the universal 
values of the transport coefficients a and k at T — ^ 0, 
and possibly their low-T corrections. We add to this 
arsenal of tools the thermoelectric coefficient, which at 
low temperatures scales linearly with temperature, with 
a non-universal slope. 
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APPENDIX A: IMPURITY SELF-ENERGY AT 
LOW ENERGIES 

The ansatz in Eq. can be used together with 

the expression for the self energy in Eq. to self- 
consistently compute the parameters a, 6, c, d, h, and 7. 
Actually, only 7 need to be computed self-consistently; 
the other parameters immediately follows. To lowest or- 



der, e , we obtain 



7 



FJi 



(Al) 



l-(Ti + cr,; J2 ' 

where Ji is the first of a scries of Fermi surface integrals 



Jn = Jnil) dp 



[A2(p^)+72] 



1/2 ' 



(A2) 



where n is an integer (only odd n appears here). In the 
next order, e^, we get 

a = rJiD-\ 

, , (A3) 

c = 27^(7,(1 -^T,:)(Ji - Js)D-\ 
where D = 2(7^7 Ji(Ji — J3) + FJ3. The next order, e^, is 
more complicated, but can be expressed in terms of the 
following functions 



^0 — Jii 



-{Ji - J3), 

7 

-(^l-^3) + ^^(J3- J5 
7 27^ 

-(^1-^3)+, „ 3 

7 \2 7 

5a^ 



3a2 
2 

1 3ab 
— 



(A4) 



{J3 - J5) 



2 J 



3 [Jb ~ -h), 



and 



mo = 1 - cr, + (Tifcg, 

2'^' h I 
mi = Kofci, 

mo 

m2 = -{kl - 2fcofc2), 

mo 

mz = -{kokz + fcife). 

mo 



(A5) 



We get 



b = D-H 2a,.h{.h - .h)a{a - 1) + — 



2(^3 - ^5)(F - 2a,7Ji) + aa{Ji - 



d = jD-^ (F - 2a,7Ji 



la' 3ab\ 5a' 

2y - — j (^3- J5)- 2:^(^5- J7) 



T^/l - a, r 3 , o 1 
n — — = [mi + 2mim2 — m3j 



/cTimo 



+ niomib + mom2(l — a) — 2crijkik2 > , (A6) 



The Fermi surface integrals can in general be computed To get explicit expressions we need a model of the or- 
through recursion (n = 1, 3, 5, ...), der parameter. For the numerics in this paper we 

J„+2 = J„-^^. (A7) 
n d'-f 
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use for simplicity A{pf) = Aq cos(2(/)j7^). in which case 

Ji = (2/Tr)VT^K[k], where k = I/a/I + (t/^o)^ 
and -fT [A:] is the complete elliptic integral of the first 



kindii In this case J3 = {2/^)^1 - k^E \k], where E [k] 
is the complete elliptic integral of the second kind, and 
the other J,i follow by recursion. 
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